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We investigate the thermal and dynamical evolution of primordial gas clouds in the 
universe after decoupling. Comparing the time-scale of dynamical evolution with that of 
fragmentation, we can estimate the typical fragmentation scale. We propose the following 
scenario of the formation process of first luminous objects consisting of large number stars. 
First, by pancake collapse of the overdensity regions in the expanding universe or collision 
between clouds in potential wells, quasi-plane shocks form. If the shock-heated temperature 
is higher than about 10 4 K, the postshock gas cools down to several hundred K by H2 line 
cooling, and the shock-compressed layer fragments into filamentary clouds. The filamentary 
cloud collapses dynamically once more and fragments into cloud cores. Finally, a primordial 
star forms in a cloud core. We show that the minimum mass of the first star is essentially 
determined by the Chandrasekhar mass. Also, we investigate the dynamical collapse of cloud 
cores by numerical simulation and show that the evolution paths of the central regions of the 
cores depend only very weakly on the total core mass. After mass accretion, a massive star 
may be formed in a core, since the estimated mass accretion rate is very large. In such a case, 
it may be possible for many massive stars form almost simultaneously. Then the clouds can 
be luminous objects. On the other hand, if the shock-heated temperature is lower, effective 
star formation is delayed significantly. 



§1. Introduction 

Today, we have a great deal of observational data concerning the early universe. 
However, we have very little information about an era referred to as the 'dark ages'. 
Information regarding a era of recombination (with redshift z of about 10 3 ) can be 
obtained by observation of cosmic microwave background radiation. After recombi- 
nation we can obtain little information until z ~ 5, when we can observe luminous 
objects such as galaxies and QSOs. On the other hand, formation stage of these 
objects should be in the 'dark ages', and thus investigating the formation process of 
galaxies using a theoretical approach is very important. 
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In the standard theory of the formation of luminous objects, the most basic 
idea is that they are made from giant gas clouds, which were once slightly over dense 
regions in the early epoch and grew to depart from the general cosmic expansion. 
After contraction of these clouds, if they fragment into many stellar-size clouds and 
many massive stars are formed on a time scale of the life of the massive stars (10 7 
yr), they can become luminous objects. To understand the way in which luminous 
objects are formed, we investigate physical processes of the clouds in various stages 
of evolution. 

To discuss the formation and evolution of galaxies, it is necessary to investi- 
gate star formation processes in protogalactic clouds, since stars are very important 
components of observed galaxies. In addition, stars play crucial roles in galactic 
activities, for example, ultraviolet radiation and supernova explosions, which are im- 
portant in some cosmological contexts. Moreover, considering the hierarchical clus- 
tering scenario, several authors have stressed the importance of stars in pregalactic 
objects for the reionization of the universe (e.g., Ref. [[])). Therefore, to study galaxy 
formation and evolution of intergalactic matter, it is necessary to investigate the star 
formation process in the early universe. 

Because gravity alone cannot be used to derive natural scales, to form ob- 
jects of various scales, such as galaxies and stars, thermal processes within gases 
are very important. In the early epoch, before the first stars were formed, bary- 
onic matter was made mostly of hydrogen and helium, and the most important 
coolants at relatively low temperature was hydrogen molecules. Without hydro- 
gen molecules, gas clouds of primordial composition are ineffective to cool below 
about 10 4 K, because of the high energy of hydrogen Lya photons. Hence the Jeans 
mass of a purely atomic gas of primordial composition is about Mj cx T 3 / 2 /? -1 / 2 ~ 
10 7 (T/10 4 K) 3 / 2 (n/l cm" 3 ) _1 / 2 M , far larger than the stellar scale. Sufficiently 
after the recombination (z iec ~ 1100), the fraction of H2 is frozen at a constant value 
of order 10 _5 ,|) and this small amount of H2 is not enough to cool the gas in order 
to lower the Jeans mass down to stellar mass within the age of the universe. Hence 
the rength scale into which the clouds fragment depends strongly on the amount of 
H2 which is formed in the clouds during their contraction. 

Because of the absence of dust grains and the inefficiency of direct reaction (H 
+ H — > H2 +7), H2 molecules are formed through the two indirect paths described 
below: 

H + e^H"+7, (1-1) 
H + H-^H 2 + e, (1-2) 



H + H+^H++ 7 , (1-3) 
H+ + H^H 2 + H+. (1-4) 

Both paths require a certain degree of ionization, and usually the first path is more 
effective than second path. Hydrogen molecules are collisionally dissociated effi- 
ciently at T ^ several thousand K. Then, in order for H2 to be made through the 
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above two reactions, there must be some degree of ionization at T 10 4 K. On 
the other hand, if in chemical equilibrium, there remain few electrons or H + ions at 
T <; 10 4 K. In such clouds, little H 2 can be formed. When z ^ 10 2 , H is radiatively 
detached by cosmic background radiation and H2 formation is inefficient. !)' W Thus, 
we consider clouds collapsing only when z & 10 2 . 

In addition, the thermal evolution of clouds is strongly affected by the initial 
shape of the clouds, or their angular momentum in the initial state. a The collapse 
of an initially oblate cloud is very different from that of spherical one. The non- 
spherical collapse in the initial phase is immediately followed by a disk-like collapse, 
which has a different time scale from a spherical collapse. Strictly speaking, the 
times needed for p to be 00 do not differ greatly from each other, but at the same 
density, which is much higher than the initial value, the time scales determined by 
the temporal shrinking rates are entirely different. For spherical collapse we have 

tdyn = —r/v r =oc p -1 / 2 (free fall time), (1"5) 

and for disk-like collapse, 

tdyn = -Z/V g (X p' 1 (z -> 0). (1-6) 

Since, after sufficient collapse, v z is almost constant and z changes as p" 1 . Then 
tdyn f° r disk-like collapse becomes much shorter than the free-fall time at the same 
density. And this time scale tdyn is also the time scale of adiabatic heating. So 
the thermal evolution of these oblate spheroids is also very different from that of 
spherical clouds. Moreover, thin disk-like clouds can fragment into smaller clouds. 
Then it is necessary to include the effect of shape change when we investigate the 
thermal evolution and the fragmentation of clouds. 

In this paper, we put together important physical processes regarding cloud 
evolution in the early universe in various evolutional stages. In §2 we investigate the 
thermal and dynamical evolution of primordial gas clouds and discuss fragmentation 
processes. In §3 we disscuss cylindrical collapse, which follows fragmentation of 
the disk-like cloud. In §4 we investigate the evolution of cloud cores and discuss 
the formation of stellar cores. Finally, we summarize our results and discuss the 
formation of first luminous objects in §5. 

§2. Contraction and fragmentation processes of primordial gas clouds 

In this section, we investigate the processes of contraction and fragmentation of 
primordial gas clouds. To discuss cloud evolution, we compare various time scales. 
And we consider that a collapsing cloud fragments when the condition 

*dyn ~ £frag> (2-1) 

is satisfied, where 




: the timescale for dynamical evolution, 



tf rag : the timescale for fragmentation. 
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2.1. Initial pressure-free collapse 

First, a cloud with a mass much larger than its own virial mass collapses almost 
as a free-fall collapse. Several authors investigated the thermal evolution of primor- 
dial clouds in pressure-free collapse with a 1-zone approximation. □''i^'i* They also 
estimated the minimum Jeans mass of a collapsing cloud. 

However, spherical clouds in pressure-free collapse are unstable against non- 
spherical perturbations & and eventually they collapse into a disk-like configura- 
tion.© The effective 'adiabatic index', r = g^J/gl , changes with the collapse con- 
figuration, and, moreover, the critical r, which corresponds to occurence of bounce, 
also changes: it is for disk collapse, 1 for cylindrical collapse, and 3/4 for spherical 
collapse^ Thus we investigated the thermal evolution of non-spherical primordial gas 
clouds. B We studied the infulence of the angular momentum and/or initial shape of 
the clouds in pressure-free collapse using a calculation of thermal evolution including 
radiative tansfer. We show that ordinary clouds bounce and are shock-heated when 
the hydrogen number density is much lower than 10 8 cm -3 , which is the density at 
which three-body reaction for H2 formation becomes important. Therefore the Jeans 
mass cannot be lowered to the stellar mass scale in this stage. 

2.2. Evolution of post schock gas 

As discussed above, ordinary pregalactic clouds collapse into disk-like configura- 
tions and are shock-heated when they bounce. Also, when subgalactic clouds collide, 
they are shock-heated. After considerable cooling, a shock-compressed layer can frag- 
ment into smaller clouds. Thus, we study the thermal evolution of the post-shock 
flow.& 

2.2.1. Important time scales 

First, we compare the time scales which are essential to the thermal evolution 
of the postshock layer in primordial gas clouds. The various relevant time scales are 



pkT 
' jim p A ' 



^cool — Cp > (2'2) 



^ lon k ion n R n e k ion n N (l-y e )' 

t iec = VjP = \ ; (2.4) 

k TQC n e n p k iec n^y e ' 

idis - 5, (2-5) 

n H 2 

**» = t||> (2-6) 
n H 2 

where t coo i) ^ioru *rec, ^dis and if or denote the time scales of cooling, ionization, recom- 
bination, H2 dissociation and formation, respectively, til njy, n«, y» (= ni/nx), k x , p, p 



*' The time scale of thermalization in the post-shock layer is shorter than the other time scales. 
For example, at 2 x 10 4 K the time scales of thermalization for plasma and neutral gas are about 
10 2 and 10 5 times shorter than the ionization time scale, respectively. 
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and m p are the number density of nucleons, the number density of the ith species, 
the fraction of ith species, the chemical reaction rate coefficient of the "X" process, 
the mass density, mean molecular weight and proton mass, respectively. A denotes 
the net cooling rate, which includes radiative and chemical cooling/heating. We 
assume that the system is optically thin for cooling photons. This assumption is 
valid for usual subgalactic clouds. For simplicity, we do not include helium. The 
existence of helium affects the mean molecular weight and the cooling rate at high 
temperature (T ~ 10 5 K). However, these are minor effects on the thermal evolution 
of primordial gas clouds at lower temperature (T 5^ 10 4 K), which we are especially 
interested in. C p denotes the heat capacity of the gas cloud. Since we are inter- 
ested in a quasi-steady post-shock flow, where i coo i is shorter than the time scale 
of the change of the shock velocity and/or pre-shock density, the post-shock layer 
is almost isobaric. and rv§i denote the effective H2 dissociation and formation 



rate. Since H2 molecu 
dissociated throgh , 



:s are formed through H~or Hg (Eqs. (IT) to (1-4)) and 
' we include only the part of reaction rate, that corresponds 
to the reaction formed from H and dissociated to H. Comparing these five time scales 
we can find the fastest process without solving detailed time-dependent differential 
equations. 

We should remark that all of these time scales are proportional to p^ 1 if the 
cooling rate and all reaction rates are proportional to p 2 . The dominant component 
of the cooling changes according to temperature. Below 10 4 K, H2 line emission 
dominates the total cooling of the cloud. In this process, the cooling rate is propor- 
tional to the number density of the hydrogen molecules at the excited levels whose 
fraction is proportional to p. LJ. Hence, the cooling time scale is proportional to p _1 
for T & 10 4 K, as in equation (2-2). The cooling process in the temperature range 
10 4 K ^ T ^ 10 5 K is dominated by the bound-bound transition of the hydrogen 
atoms, and the number of excited hydrogen atoms is determined by collisions with 
other atoms. Therefore, the cooling rate is proportional to the square of the total 
density. For higher temperatures (T 10 5 K), the free- free emission by the collision 
between ions and electrons dominates the energy radiated away from the cloud. In 
any case^JJie cooling rate is proportional to p 2 , and the cooling time is proportional 
to bll The chemical reaction rates are proportional to p 2 if the reactions are 
dominated by collisional processes. 

Photoionization and H2 photodissociation processes caused by the UV photons 
emitted from the post-shock hot region could affect the ionization rates and the H2 
dissociation rates in post-shock regions. These radiative reactions are proportional 
to p 2 , because the emission rate of UV photons in the post-shock is pronortional to 
p 2 . We investigate the effects of UV photons from postshock regions. i3 However, 
as a Oth-order approximation, we neglect the photoionization and the H2 photodis- 
sociation processes in this paper. 



*' This treatment is limited to low density (n 5s 10 4 cm -3 ) The cooling rate is not proportional 
to p 2 for higher density (higher than the critical density). 

Compton cooling is effective at high redshift and high temperature. In this case, the cooling 
rate is not proportional to p 2 , but this does not change the following results significantly. 
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2.2.2. Fraction of H 2 

Here we consider the estimation of yn 2 m the post-shock flow for the case in which 
the post-shock temperature is higher than about 10 4 K. If yn 2 can be determined for 
given T and y e , the cooling time is always determined by just T and y e . As a result, 
the thermal evolution of the system is completely determined in the y e -T plane. 

Now, we make some remarks regarding important properties of tdis and tf or - 
First, the time scale of dissociation is independent of 2/h 2) the fraction of hydrogen 
molecules. On the other hand, tf or is proportional to yn 2 - Then, the time scales 
given in terms of their equilibrium values are 



*dis - *dj s > ( 2 ' 7 ) 
, _ ,eq VU 2 _ ,eq 2/H 2 (r> n\ 

''" r "l/uM/.-T) S&(v.,T)' 



where the suffix eq denotes the value at which the hydrogen molecules are in chemical 
equilibrium for a given electron abundance and temperature. Thus, idis < ^reo *cooi 
is a sufficient condition yn 2 to be the equilibrium value for given y e and T. 



t e 1 ,,t >t.. E 
cool' rec dis H 
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Fig. 1. The y e -T plane is divided into two regions, one in which chemical equilibrium is achieved 
for a fraction of hydrogen molecules and one in which this is not the case. The solid boundary 
of the shaded region is determined by the condition t coo \ = t&s. The long dashed line denotes 
the boundary at which t ICC = tdis is satisfied. In the shaded region, the chemical equilibrium of 
hydrogen molecules breaks down, and the fraction j/h 2 freezes as the temperature and/or ion- 
ization degree drops. The dotted lines labeled '100 km/s' and '20 km/s' denote the evolutionary 
path of j/ e . 

We compare tdis with t rec and t COQ \ in the y e -T plane (Fig. 1). At higher 
temperature, tdis is shorter than t rec and t coo \, and then yn 2 = y^Aye,T). After 



*' Note that Fig. 1 is corresponds to Fig. 2 of Susa et al.,LP but calculated using new reaction 
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cooling, in the shaded region, the chemical equilibrium of hydrogen molecules breaks 
down, and yn 2 almost freezes as the temperature and/or ionization degree drops. 
Thus, we can roughly determine yn 2 f° r given y e and T in the cooling layers of 
post-shock flow. 

2.2.3. Ionization and cooling 

Here we investigate the case t; on < t rec , which corresponds to higher temperature 
(T <; 10 4 K), and compare the time scale of ionization with that of cooling. The 
line along which ti on /t coo \ is equal to unity is indicated by the solid line in Fig. 2. If 
the initial conditions are given below the line, the system is ionized before it cools 
down, because tionAcool < 1- Otherwise, the system cools before it is ionized. The 
expected evolutionary path is also expressed by the arrows in Fig. 2. 

2.2.4. Recombination and cooling 

Next we investigate the case t Tec < ti on , which corresponds to lower temperature 
and examine the balance between the recombination process and the cooling process. 
The recombination process is important for T 5s 10 4 K . At T 5s 10 4 K, the 
line cooling of the hydrogen molecules dominates the atomic hydrogen line cooling. 
Therefore the ratio t rec /t coo i depends not only on the ionization degree and the 
temperature but also on the fraction of hydrogen molecules. In order to draw the 
line irecAcool = 1 in the y e -T plane, we need information on the fraction of hydrogen 
molecules. The fraction of H2 is determined by the chemical equilibrium above 
~ 10 4 K. Then, it is approximately frozen below ~ 6x 10 3 K through an evolutionary 
course with cooling. The line along which t Tec /t coo \ = 1 is satisfied is drawn in Fig. 2. 
In the region above the line, the system recombines before it cools down, and below 
the line, the system cools down before the recombination proceeds. The expected 
evolutionary path is also indicated in Fig. 2 by the arrows. 

2.2.5. Shock diagram 

The ratio of any two time scales is independent of the total density, and they are 
determined by the temperature and the chemical compositions (y e and yn 2 ). yn 2 1S 
almost always determined by y e and T. Therefore, the ratio of any two time scales 
just depends on T and y e . In the following sections, we compare the time scales 
individually. 

The main results are summarized in Fig. 2. The evolution of a shock heated 
system is basically explained in the y e -T plane, which we call "shock diagram". 
The evolutionary path of the post-shock layer is obtained by tracing the directions 
of the arrows on the shock diagram. Post-shock gas, which is heated up to high 
temperature, appears in region (a). In region (a), the shortest time scale is ti on , 
and the system evolves upward and enters region (b). In region (b), i coo i is the 
shortest, and therefore the system evolves leftward here and comes into region (c). 
In region (c), i rcc is the shortest. Then the system evolves downward and reenters 
the region (b), where it evolves leftward again and crosses the boundary of regions 



rates and cooling function. B 1 There exist some quantitative diferences, however. In particular, the 
line tree = £dis changes somewhat. This is mainly because the dissociative recombination rate of 
(H2 + + e — > 2H) is reduced,™ and hence the effective H2 dissociation rate is reduced. 
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Fig. 2. The y e -T plane is divided into several regions in which the ratio of the time scales are 
different from each other. The arrows on the plane denote the directions of the evolution of the 
post-shock layer. The dot-short dashed line denotes the curve on which the equality t TCC — ti on 
holds. The region (a) in which the condition ij on < t coo \ holds is bounded by the solid line. In 
this region, the system evolves upward in the y e -T diagram. In the region (b), the conditions 
tcooi < tion, ircc always holds. Therefore the temperature drops down before the system is 
further ionized or recombines. In the region (c) we have t ICC < t coo \, and the system evolves 
downward. The region (c) is bounded by the dotted line, short dashed line and the long dashed 
line. These three lines denote the boundary on which the condition i coo i = t lac holds. The 
long dashed line is calculated by assuming that the cooling rate is dominated by hydrogen 
line cooling. The dotted line is calculated assuming the chemical equilibrium of the hydrogen 
molecules, and the short dashed line is calculated with ya 2 = 10 -3 , which corresponds to the 
frozen value of the hydrogen molecules for a shock velocity higher than ~ 40 km/s. 

(b) and (c). Finally, the system evolves nearly along the short dashed curve which is 
the boundary of regions (b) and (c). For post-shock gas which is not heated above 
2 x 10 4 K, the thermal evolution is somewhat different from the previous case. In 
this case, when the system enters region (b) from region (a), the ionization degree is 
so low that the system does not enter region (c) at T ~ 10 4 K. Then the ionization 
degree and the fraction of hydrogen molecules do not "forget" the initial conditions 
of the chemical composition. 

2.2.6. Convergence of y e and yn 2 111 almost steady postshock flow 

It is known that in steady post shock flow, both y e and yn 2 become the same 

value at the same temperature. This is surprising because ionization equilibrium is 

not yet achieved at all. Now, the reason that covergence of y e and yu 2 in the steady 

post-shock flow occurs can be understood using Figs. 1 and 2. 

At T ~ 10 4 K, the line i rec = t coo \ becomes nearly vertical, since the dominant 

cooling process is the hydrogen atomic line cooling, which decreases by many orders 
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below T ~ 3 x 10 4 K. If (1) the initial shock velocity is larger than ~ 40 km/s or 
(2) the initial ionization degree is larger than ~ 5x 10 -2 and the initial temperature 
is larger than T ~ 10 4 K, then the evolutionary path should hit this nearly vertical 
line (the long dashed line in Fig. 2). Because the line t rec = t coo \ is nearly vertical 
at T ~ 10 4 K and y e ^ 5 x 10~ 2 , y e drops to several xl0~ 2 without changing 
the temperature, and "forgets" the initial conditions, though y e is not in chemical 
equilibrium. The value to which y e converges (~ 5 x 10~ 2 ) is essentially determined 
by the equation t iec (y e ,T) = t^ ol (y e ,T) at T ~ 8000 K. Here, the temperature 
8000 K is obtained by the equality t rec /t^°^ = 1, where is the cooling time 
estimated by the hydrogen line cooling. The equation t iec /t^°^ = 1 depends both 
on y e and on T in general, but it depends only on T for y e <C 1. 

The reason for the convergence of H2 also can be understood. According to 
Fig. 1, H2 is still in chemical equilibrium just below T = 10 4 K, at which the 
convergence of y e takes place. As a result, the systems that satisfy the previous 
condition (1) or (2) experience the same state (same y e , same yn 2 ) just below T = 
10 4 K. Therefore yn 2 also converges for different initial conditions. 

2.2.7. Comparizon with numerical calculations 

1 

0.1 
0.01 
0.001 

V 

0.0001 



10- 5 
10- 6 
10- 7 

100 1000 10 4 10 5 10 6 
T [ K ] 

Fig. 3. The evolution of the post-shock layer in the y e -T plane. The solid lines are the boundaries 
beyond which the direction of evolution is changed, and each line has the same meaning as in 
Fig. 2. The short dashed lines are the numerical results. Each suffix denotes the initial shock 
velocity. 

In this section we compare the results obtained in the previous section with 
numerical calculations. Figure 3 displays the evolution of y e as a function of the 
temperature behind the shock front in steady flow. It is obvious that the expected 
evolutionary path in Fig. 2 agrees with the numerically calculated path. We also 
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find that there exists a critical value of the initial shock velocity above which the 
evolution of y e below 10 4 K converges. The critical value of the shock velocity is 
~ 40 km/s. This condition for the initial shock velocity is equivalent to the condition 
(1) in the previous subsection. 

Because of the delay of the recombination process compared to the line cooling 
due to the hydrogen molecules, fairly high degree of ionization (y e ~ 10~ L5 ) is possible 
below 8000K. Feeded the electrons, H2 also maintains a high abundance {yu 2 ~ 
10 -2 ) for T Si several x 10 3 K, then the cloud temperature drops to T ~ 100 K 
quickly due to the H2 line cooling. In fact, the predicted evolutionary paths agree 
well with the numerically integrated paths for steady flow. 

2.2.8. Fragmentation of the post-shock layer 

To estimate the fragmentation condition of the post-shock layer, we can consider 
t C ooi as tdyn and ts as if rag , where ts is the free-fall time scale, which coresponds 
to the growth time scale of gravitational instability. HHP With the^standard CDM 
cosmol og ical model, we can estimate the fragmentation epoch. EJ) As shown in 
Fig. 4,li 3 l primordial gas clouds, which collapse in the early universe, fragment after 
they are clooded down to about 100 K if M J> 1O 8 M . In the case that M <, 1O 8 M , 
the post-shock temperature is lower than 10 4 K, and H2 is formed only slightly. So, 
from the initial state in which t coo i > tg, the cloud evolves by gravitational force 
without cooling. Then the cloud will be virialized and cannot fragment, because it 
cannot be thin enough for fragmentation. 

§3. Collapse and fagmentaion of filamentary clouds 

A cooled sheet fragments more easily into filaments!!!''© than into spherical 
clouds. As the virial temperature of a filament is essentially constant, its evolution 
is much different from that of a spherical cloud, whose virial temperature increases 
as the cloud contracts. In this section, we investigate the evolution of filamentary 
clouds of primordial gas and estimate the fragmentation epoch. 

3.1. Basic equations 
We investigate 



;he dynamical and thermal evolution with a 1-zone approximation 
for the simplicity. tZP Then, basic equations are as follows. For an isothermal cylinder 
there is a critical line density M C (T) given by 

A « T » - ^ < 3 "> 

where T, [/,, and mn represent the isothermal tenrnerature, the mean molecular weight 
and the mass of hydrogen atoms, respectively.^ An isothermal cylinder with line 
density M larger than M C (T) has no equilibrium and collapses, while an isothermal 
cylinder with M less than M C (T) expands if the external pressure does not exist. 

When the filamentary cloud is formed as a result of the fragmentation of the 
sheet, the gas pressure is comparable to the gravitational force. Hence the effect 
of the gas pressure must be included in the equation of motion of the cylinder like 
cloud. To know the dynamical evolution of such a cylinder we begin with the virial 
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Fig. 4. Estimated fragmentation epoch. Pre-shock density-shock velocity plane, which is divided 
into three regions. Region A is where fragmentation occurs at low temperature ~ 100 K, and 
region B is where it occurs at several thousand kelvin. In region C the post-shock temperature 
cannot rise above several thousand kelvin. Typical values of the shock velocity and the pre- 
shock density at the collapse in the standard CDM model are overplotted (dashed line). Numbers 
plotted above this line indicate the mass scales in unit of solar mass. The dot-dashed line divides 
the plane according to whether fragmentation occurs within 1 Gyr. The time elapsed since the 
point at which the shock front cross until fragmentation is calculated with z = 5. Most of the 
galactic sized clouds are cooled within 1 Gyr. Adapted from Yamada and Nishi (1998). 



equation for a cylinder, 



where 



1 d 2 I 

=2T + 277-GM 2 , (3-2) 



7 = / pr 2 dV, T = [ \pu 2 dV, 11= [ pdV, (3-3) 
Jv Jv 2 Jv 

and the integration is effected over the volume per unit length of the cylinder, V . 
Here we assume that the cylinder is uniform and has constant density p, scale radius 
R = ^M/irp, and temperature T. Then /, T and 77 introduced above become 

' W, T = -m(^)\ n = ^M. (3-4) 



2 4 V dt J ' prriH 

Substituting Eq. (3-4) in Eq. (3-2) we obtain 

d 2 R 2G 



= - ( M -M C (T)). (3-5) 



dt 2 R 
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We describe the evolution of the cylinder by this equation. 

The time evolution of the cloud temperature is described by the energy equation 
given by 

-77 = ~P~TL ^rad ~ A:hem , (3-6) 

at at p 

where e, yl ra d and ^lchem are the thermal energy per unit mass, the cooling by the 
radiation, and the cooling/heating by chemical reactions, respectively. 

Assuming that the collapsing cylinder fragments when the condition (2-1) is 



satisfied, we estimate the epoch of fragmentatioin and the mass of fragments. The 
estimated mass scales correspond to the stellar mass scale (1M Q ~ WOMq). 

3.2. Analytical estimate of the minimum fragment mass 

Now, we estimate the minimum fragment mass analytically. As the cylinder 
collapses isothermally as a function of time in the later stage, we can set de/dt ~ 
in Eq. ( |3-6| ) and from this we obtain t^yn ~ (7 — l)*cooh where 7 = 7/5 for a diatomic 
molecular gas. Since the cloud cools by the line emissions of hydrogen molecules and 
is optically thick to these line emissions, t coo \ can be estimated by 



1 M 



c ° o1 2vRcrT^a c ' { ] 

where a = 2iT 5 k B /15h s c 2 is the Stefan-Boltzmann constant, and a c is the effective 
number of line emissions. Since the optical depth is ~ 100 at most, the wing of 
the line profile may not be important, and the line broadening is determined by 
Doppler broadening as Avjv = uh 2 / c = V^bT '/muc 2 . The temperature of the 
cylinder is estimated by the virial temperature as 

k B T = ^fimnGM. (3-8) 

From the above equations we can estimate the fragment mass as 

M,, V ~2.RM~JI-L^, (3-9) 



where mpi = y/hc/G is the Planck mass. Equation (3-9) implies that -/Wf rag is 
essentially equal to the Chandrasekhar mass (~ mpi/rrip, where m p is proton mass). 

To this point we have not considered other heating processes, such as shocks and 
turbulence, which likely occur in a collapsing cloud. All these processes tend to halt 
the collapse, so that the epoch of fragmentation becomes earlier, and, as a result, 
the mass of fragments increases. Thus Eq. (3-9) gives a lower limit on the mass 
of the fragment which is formed from a primordial gas cloud. Our estimate shows 
that even a primordial gas cloud which cools most effectively by hydrogen molecules 
cannot fragment into a smaller mass than the expressed in Eq. (3-9). 

For smaller line densities the gas pressure force halts the collar>se before three- 
body reactions completely convert atomic hydrogen into molecules. B> In these cases, 
Eq. (3-9) is not applicable, because clouds are optically thin to line emissions if the 
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molecular fraction is smaller than about 10 , and chemical heating by H2 formation 
must be considered when hydrogen molecules are formed by three-body reactions. 
These clouds fragment at lower density, in other words at larger scale radius. Since 
Mf rag is proportional to the scale radius at fragmentation, this leads to the larger 
fragment mass. 

In reality, numerical simulations suggest that primordial cylindrical clouds frag- 
ment at lower densit y, tynd the estimated masses of fragmets are larger than the 
Chandrasekhar mass.M'Ej) 



§4. Evolution of cloud cores 



Now we study the evolution process of cloud cores. After they collapse into 
stellar cores, true stars can form. 

4.1. Quasi-static contraction 

The manner in which primordial clouds of stellar mass, which are fragments 
of filamentary clouds, evolve after the fragmentation depends on their efficiency of 
cooling, i.e., the ratio of the free-fall time scale tg to the cooling time scale t coo \. 
If ig < t C ooi (i-e., cooling is not effective), while the cloud cools it has enough time 
to adjust itself to a new hydrostatic configuration. Then it contracts, maintaining 
a nearly hydrostatic equilibrium (i.e. Kelvin-Helmholtz contraction). On the other 
hand, in the case tg > t coo \, the cloud collapses dynamically in its free-fall time scale 
after substantial cooling, and no hydrostatic equilibrium is established. 

At the time of the fragmentation of filamentary clouds >-the two time scales of 
fragments, i coo i and tg, are of the same order of magnitude. E3) Thus we investigate 
quasi-static contraction of cloud cores as a first step.EiP 

4.1.1. Radiative processes in primordial molecular cloud cores 

We consider spherically symmetric gas clouds with stellar mass (~ 1M & ) ('molec- 
ular cores') composed only of H2 molecules, because the hydrogen atoms are con- 
verted into molecular form before the fragmentation for these low mass fragments. £3 
We can neglect helium, since it is thermally inert at these low temperatures. As men- 
tioned above, these clouds lose their thermal energy by radiative cooling via the H2 
rotovibrational lines, and presumably contract into protostars. 

In this section, we describe our calculation scheme for the luminosity and cooling 
rate of a spherically symmetric cloud. The specific intensity I u (ergs sec -1 ran -2 sr _1 
Hz^ 1 ) along a ray is calculated by solving the radiative transfer equation,^ 

*k = - avIv+jv . (4.1) 

ds 

Here s is the displacement along the ray, and a v and j u are the absorption and 
emission coefficients, respectively. These coefficients can be written using Einstein 
A- and .B-coefficients: 

ju = j ^n 2 A 21 ^){v) ) (4-2) 
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ot v = ^-(f>{v){niBi2 - n 2 B 2 i). (4-3) 

47T 

Here n\ and n 2 are the number densities of the molecules in the lower and upper levels 
of the transition, respectively, and 4>(y) is the line profile function. In our calculation, 
the line center optical depth of the lines across a cloud r c is no more than about 
several hundred. Since the Lorentz wings of a line become more important than its 
Doppler core only when r c is as large as 10 3 , we only consider line broadening owing 
to the thermal Doppler effect. 

Note that the critical number density n cr , where radiative and collisional de- 
exicitation rates become equal for molecular hydrogen, is ~ 10 4 cm" 3 . Thus, for 
a typical number density of fragments (~ 10 10 cm~ 3 » n cr ), almost all excited 
molecules would be deexcited by collision with other molecules. Therefore, colli- 
sional local thermodynamic equilibrium (LTE) is established, and the scattering of 
photons can be neglected. These properties make the transfer problem of H2 lines 
fairly tractable. 

4.1.2. Evolution of fragments 

As initial conditions, we assume polytropic gas spheres in hydrostatic equilib- 
rium; namely the density/temperature distribution is represented by the Emden 
function of polytropic index N. Clouds are cut off at the radius r s , where the den- 
sity falls off by a factor of 10~ 3 from the central value. Therefore, the parameters 
characterizing initial conditions are the effective polytropic index N, the total mass 
of the cloud M, and the number density at the half mass radius. For a given 
initial configuration, we can obtain the specific entropy distribution s(m,t = 0). 

We calculate quasi-static contraction as follows. First we calculate the cool- 
ing rate A(m), and advance the specific entropy distribution s(m,t) using the heat 
equation 

§ = ~A{m). (4-4) 

With the new entropy distribution, we find the new equilibrium configuration using 
the equations of hydrostatic equilibrium. 

For given states of clouds, we compare two time scales, the free-fall time 



and the time scale of quasi-static contraction, 

t qsc (m) ee pi . (4-6) 

\<JZ / m , quasi-static 

Here, p(m) is the mean density at the mass coordinate m, and each quantity is 
evaluated at the coordinate m. From the virial theorem, we know that t qsc is of the 
same order of magnitude as t coo \. 

Figure 5 displays the contours of the ratio tg/t qsc at the center (m = 0) of 
clouds with (a) iV = 2.5 (adiabatic stratification) and (b) iV = 8. Regions around 
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Fig. 5. Contours of the ratio tff/i qsc at the centers of clouds of polytropic index (a) N = 2.5 and 
(b) N = 8. The contour spacing is logarithmic with an increment of 0.5. The upper-right 
regions correspond to temperatures T > 2000 K, where hydrogen molecules begin to dissociate. 
Adapted from Omukai et al. (1998) 



the curves tg/t qsc = 1 are expected to be initial states of fragments. Although Fig. 5 
displays the values at the center, this ratio does not change significantly within a 
cloud, except near the surface. The upper-right regions correspond to higher central 
temperature > 2000 K, where the molecular hydrogen begins to dissociate. From 
the two panels of Fig. 5, we can see that the ratio tg/t qsc depends only slightly on 
the effective polytropic indices N, rather than on masses and densities. Note that 
for clouds of the same mass, the ratio tg/t qsc is greater for the denser cloud. This 
implies that as a cloud contracts, the ratio tg/t qsc becomes larger. In particular, 
for fragments of filamentary gas clouds, for which iff ~ t qsc initially, t qsc becomes 
shorter than tg as they contract; i.e., such clouds collapse dynamically. Note that 
when t qsc < tg, clouds do not contract quasi-statically, and t qsc does not possess the 
meaning of the collapse time scale. In this case, t qsc merely measures the time scale 
of cooling. 

Figure 6 exhibits the change of time scales and their ratio for the quasi-static 
contraction of a cloud which is initially at the asterisk position of Fig. 5. We can 
see that the cloud indeed collapses dynamically after it contracts quasi-statically to 
some extent. 

4.1.3. Line profile and cooling rate 

We have shown that the stellar mass molecular cores collapse dynamically, al- 
though some lines are optically thick at the line center. This fact indicates that the 
efficiency of cooling by H2 lines does not fall significantly even in such a situation, in 
contrast to the case of continuous radiation. The reason for this can be interpreted 
as follows. 

Figure 7 displays the change of a processed profile of a single optically thick 
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Fig. 6. Evolution of the free-fall time ts and dynamical time for quasi-static contraction i qsc and 
their ratio at the center for a cloud which contracts quasi-statically at the beginning (asterisk of 
Fig. 5). The time scales tg and t qBC and the time ellapsed from the begining t are all normalized 
by iqsc.o, the initial value of t qsc . This figure indeed shows that the cloud begins to collapse 
dynamically after ~ 1.5 i qsc ,o- Adapted from Omukai et al. (1998) 




-4 -2 2 4 



{v-v Q )/ Ai/ D 

Fig. 7. Processed line profiles of an optically thick line along a ray passing through the center of a 
cloud at radii containing 0% (center, solid line), 1.4% (dotted line), 9.9% (short dashed line), 
and 100% (surface, long dashed line) of the total mass. The optical depth at the line-center 
frequency r c = 38, 46, 54 and 76, respectively, measured from the surface of incidence. The 
line corresponds to the transition from (v = 1, j = 3) to (v — 0, j = 1). The parameters of 
the cloud are the same as in Fig. 2. The abscissa represents the frequency from the line center 
normalized by the Doppler width at the center of the cloud Avu,c- Adapted form Omukai et al. 
(1998) 
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line (optical depth at line center across the cloud is 76) along a ray passing through 
the center of the cloud. Since the cloud is optically thick at line center, the specific 
intensity at line center is saturated with the blackbody value determined from the 
local temperature. On the other hand, the optical depth at the wing is smaller than 
unity. The final processed profile at the cloud surface is double peaked, and the 
peak frequencies correspond to frequencies where the optical depth across the core 
TV, core > 1 (he., the radiation is saturated with the value of the core temperature) 
and that of envelope r^ env < 1 (i.e., the radiation can be larger than the value 
of the envelope temperature). The height (i.e., the normalization of intensity) and 
width of these peaks are of the same order of magnitude as the blackbody value 
determined from the temperature near the center, not around the surface of the 
cloud. The reduction factors from the central values are smaller than about 10. 
Owing to the effect described above, which is a distinctive feature of line cooling, the 
temperature around the cloud center can be 'seen' from outside the cloud, even when 
lines are optically thick at line center and the surface temperature is relatively low. 
Accordingly, its cooling proceeds efficiently, as long as there are enough frequency 
ranges where T V)COie > 1 and r^ env < 1. Therefore, the decrease in the efficiency of 
cooling caused by the low surface temperature is not severe in the line cooling case, 
in contrast to the case of cooling by continuum. 

4.2. Dynamical collapse of primordial protostellar clouds 

As discussed above, although the cooling time of a cloud is longer than the its 
free-fall time initially, the cooling time at the center becomes shorter than the free-fall 
time, and the cloud begins to dynamical collapse after some quasi-static contraction. 
In this section we investigate the dynamical collapse of primordial protostellar clouds 
and the formation of stellar cores. c3 

4.2.1. Method of calculation 

We assume spherical symmetry for the simplicity. Then the Lagrangian equa- 
tions describing dynamical collapse are as follows. 
Equation of continuity: 



equatioin of motion: 



equation of energy: 



equation of state: 



Dv _ 2 dp Gm 



" -3", (4-8) 



Dt dm 



r 



De D fl\ A . . 

Dt ~ -»Di U ~ ~ P ' (4 ' 9) 



P = bad - l)f*. (4-10) 



In the above expressions, m, p, v, p, e, A and "y a d are the mass within a radius r, the 
density, velocity, pressure, thermal energy per unit mass, net cooling rate per unit 
volume and the adiabatic coefficient, respectively. 
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The net cooling rate A consists of two parts, the radiative part yl ra d an d the 
chemical part A:hem- ^rad can be written as 

A Iad (m) = p^-, (4-11) 
am 

where L(m) is the luminosity and is obtained by solving radiative transfer equa- 
tions. We consider H2 lines and continuum components which consist of H2 collision- 
induced absorption (CIA), H~ bound- free absorption, Lyman continuum absorption, 
etc. as sources of opacity. /t c hem is given by 

, f \ de chcm (m) 

Ahcm(m) = p — , (4- 12) 

where e c h em is the chemical binding energy per unit mass. 
4.2.2. Dynamical evolution of protostellar clouds 

The evolutionary sequences of the radial profiles of number density, temperature, 
velocity and H2 concentration distributions are illustrated in Fig. 8. 

The collapse proceeds almost self-similaly and is analogous to the Larson-Penston 
similarity solution© until the central region reashes the stellar density. The dashed 
line in the figure indicates the slope of density gradient of —2. The slope in this 
case is slightly steeper and its value is aboud —2.2, which correspond to the effective 
adiabatic index is 1.1. 

At higher density, almost all hydrogens are converted to H2 by efficient three- 
body processe. But the total mass of central molecular part is about 1 M & which 
hardly depends on total cloud mass. When the central number density reaches about 
3 x 10 13 cm -3 and the central temperature reaches about 1600 K, the hydrogen 
molecules begin to dissociate gradually. But at almost the same time, H2 CIA 
continuum cooling becomes efficient. This continuum cooling is so strong to stop 
and even to reverse dissociation. After a small central part of the cloud becomes 
opaque to CIA continuum (3 xlO 16 cm' 3 , 2000 K), the effective dissocication begins. 

The evolution of the temperature and the effective 'adiabatic index' r = g^/^ 
at the center for various total mass cases are illustrated in Fig. 9. We can see the 
central value r is about 1.1 almost throughout the collapse. After some contraction, 
evolution paths of the central regions are hardly affected by the total mass. 

When most of the H2 molecules are dissociated, r rises above the critical value 
4/3. After the central part of the cloud contracts almost adiabatically to some 
extent, a hydrostatic core with very small mass (~5x 1O~ 3 M0) forms at the center. 
At that time, the central number density and temperature are ~ 10 22 cm -3 and 
~ 3 x 10 4 K. These values are almost the same as those of present-day stellar cores 
at the formation epoch. 

But before the formation of the stellar core, no transit core is formed, in contrast 
to the case of present-day star formation. And the estimated mass accreation rate 
in the case of primordial star formation is much different from that of present-day 
star formation (~ 1O~ 5 M0 yr _1 ), as shown below. Befor stellar core formation, the 
radial profiles of the density are well fitted by the Larson-Penston-type similarity 
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Fig. 8. The radial profiles of the physical quantities, (a) The number density (n = p/mn, where mn 
is the mass of the hydrogen atom), (b) temperature, (c) velocity, and (d) the H2 concentration 
are given as functions of the radial distance. The numbers in the figure denote the evolutionary 
sequences. '0' corresponds to the initial time, and '1' to 5.7 x 10 5 yrs after starting calculation. 
Here, three-body processes are active in the central region and temperature inversion occurs 
because of efficient H2 line cooling. '2'corresponds to 8.5 x 10 3 yrs after state 1: The cloud 
becomes optically thick to some H2 lines. '3' corresponds to 2.8 x 10 2 yrs after state 2: The 
central region becomes fully molecular. '4' corresponds to 12 yrs after state 3: The central 
region become optically thick to the H2 CIA continuum. '5' corresponds to 0.32 yrs after state 
4: H2 dissociation becomes efficient. '6' corresponds to 2.4 x 1CF 2 yrs after state 5: Shortly 
after the core formation. '7' corresponds to 4.1 x 10~ 2 yrs after state 6. Adapted from Omukai 
and Nishi (1998) 



solution with 7 = 1.09. Assuming that the evolution of an envelope after stellar core 
formation is described by the same similarity solution, we can estimate the evolution 
of the core mass. With the similarity solution after the central density becoming 
infinite, we have for the central stellar massif 

M* = O.11M 1 ( 4 -!3) 
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Fig. 9. The temperature and effective 'adiabatic index' r at the central regions are plotted against 
the central number density n c [cm ] for clouds with various masses, (a) The central tempera- 
ture T c . (b) The effective 'adiabatic index' r c in the central regions. In both panels, the total 
masses are 1O 2 M (solid lines), 10M Q (dotted lines) and 1M (dashued lines). Adapted from 
Omukai and Nishi (1998) 

and we can estimate the mass accretion rate as 

/ + \-0.27 

M* = 8.3 x KT 2 M yr^ 1 J , (4-14) 

= 3.7 x 1O~ 2 M y^ 1 

The mass accretion rate is very large and diminishes with time. The main reason 
for this large mass accretion rate is that the temperature of primordial gas clouds 
(~ 10 3 K) is higher (which is the result of weak cooling) than that of the present-day 
molecular clouds (~ 10 K). The final stellar mass may be large because of the large 
mass accretion rate. 

§5. Summary and discussion 

We investigated the thermal and the dynamical evolution of primordial gas 
clouds in the universe after decoupling and discussed the fragmentation process. 
It is necessary to study the evolution of the clouds considering the change of cloud 
shapes. Comparing the time scale of the dynamical evolution with that of fragmen- 
tation, we estimated the typical fragmentation scale. 
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5.1. A formation scenario of first luminous objects 



We propose a formation scenario for first luminous objects (Fig. 1C) 



CDM 
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baryon 
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dark halo 
shock 



cloud core 




f ragmentat i on 

i nto 
cy I i ndr i ca I 




Fig. 10. The formation scenario for first luminous objects. 

(1) By pancake collapse of the over-density regions or collision between clouds in 
potential wells, quasi-plane shocks form. 

(2) If the shock-heated temperature is higher than about 10 4 K, the post-shock 
gas cools by H2 line cooling, and the shock-compressed layer fragments into 
cylindrical clouds when td yn ~ ifrag (T is several hundred K). 

(3) The cylindrical cloud collapses dynamically once more and fragments into cloud 
cores when i^yn ~ ^frag- 
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(4) A primordial star forms in a cloud core. 

(5) In the case that the shock-heated temperature is higher than about 10 4 K, many 
fragments can be formed almost simultaneously. Thus, many massive stars may 
be formed within the time scale of the life of massive stars. 

(6) Luminous objects are formed. 

The condition (5) is achieved when a massive cloud (<; 1O 8 M0) collapses or 
subgalactic clouds collide with each other in the potential of massive objects, which 
has the virial velocity larger than about 40 km/s. In any case, luminous objects can 
be formed after massive objects (M ^ 10 8 M Q ) collapse. 

In the case that the shock-heated temperature is lower than about 10 4 K, ioniza- 
tion can not occur and H2 molecules are formed only with cosmological relic electrons. 
Hence yn 2 is small and cooling by H2 is weak, and the cloud cannot fragment into 
many small clouds. In this case, one or a few massive stars may be formed around 
the cloud center. After the first generation of massive stars are formed, the cooling 
of the remaining part of the cloud becomes inefficient, since hydrogen molecules are 
dissociated by UV photons radiated by massive stars. Thus, the star formation rate 
is strongly self-regulated.© Moreover, the binding energy of such a small mass cloud 
is small, and the cloud is blown off or disturbed severely by a super nova explosion, 
and as a result, the formation of the next generation of stars is greatly delayed. Then 
it is unlikely for isolated small mass clouds to become luminous objects. 
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